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The momentum-shell renormalization group (RG) is used to study the condensation of interacting 
Bose gases without and with disorder. First of all, for the homogeneous disorder-free Bose gas the 
interaction-induced shifts in the critical temperature and chemical potential are determined up to 
second order in the scattering length. The approach does not make use of dimensional reduction 
and is thus independent of previous derivations. Secondly, the RG is used together with the replica 
method to study the interacting Bose gas with delta-correlated disorder. The flow equations are 
derived and found to reduce, in the high-temperature limit, to the RG equations of the classical 
Landau- Ginzburg model with random-exchange defects. The random fixed point is used to calculate 
the condensation temperature under the combined influence of particle interactions and disorder. 



PACS numbers: 03.75.Hh,05.30.Jp,64.60.Ak 

The realization of Bose-Einstein condensation in ul- 
tracold atomic gases has renewed the interest in the 
quantum-statistical properties of interacting Bose gases 
around the phase transition. One particular focus of the- 
oretical work in this area concerns the calculation of the 
critical temperature T c of the homogeneous Bose gas as a 
function of the scattering length a. Due to its nonpertur- 
bative character, this apparently simple problem turned 
out to present a considerable theoretical challenge. For a 
long time, the leading-order functional dependence of the 
shift in T c on the scattering length remained highly con- 
troversial. It was even unclear whether the interactions 
lead to an upward or downward shift in the critical tem- 
perature p|. Only recently, a generally accepted result 
has emerged 0, tJ El • It was shown that the critical tem- 
perature of an interacting Bose gas with spatial density 
n obeys the relation 






1 + c^n 1 ' 3 + [c' 2 Han 1 / 3 ) + c'^n 2 ' 3 
+0{a 3 n), 



(1) 



where Xj? denotes the critical temperature of an ideal 
Bose gas with the same density. To calculate the coeffi- 
cients Ci, one can make use of dimensional reduction 1, 5] 
and perturbatively match the full quantum-mechanical 
problem to a classical field theory with a lower- (i.e., 
three-) dimensional action |j]. The coefficient c' 2 can then 
be determined by perturbation theory alone, whereas the 
evaluation of c\ and c 2 ' requires nonperturbative input 
from the classical theory. This input was calculated in 
Refs. [fluHS! with so far unsurpassed accuracy by means 
of Monte Carlo lattice simulations. 

Even after the qualitative and quantitative behavior of 
the T c shift had been solidly established US SHU, a 
large number of papers continued to appear that pre- 
sented further calculations of critical properties of in- 
teracting Bose gases. The main purpose of these works 
was not to correct or improve the previous results, but 
to introduce, test, and refine alternative approaches. 
The a pplie d methods comprise, e.g., linear S expansion 
IS 113 Hi El El, variational perturbation the- 




and the exact renormalization group 



The present paper follows this research line of inves- 
tigating alternative approaches to studying the interact- 
ing Bose gas at criticality. One of its main motivations 
is the observation that most of the methods that have 
been used so far, e.g., the perturbative matching proce- 
dure of Ref. |4| for finding the second-order shifts as well 
as the nonperturbative calculations within the classical 
theory, are relatively involved. This raises the question 
of whether there exist alternative approaches that, on 
the one hand, are more straightforward to apply, but, 
on the other hand, still yield sufficiently accurate quan- 
titative results and are also readily applicable to more 
complicated problems. In the present paper it is shown 
that the one-loop momentum-shell renormalization group 
(RG) method J2l| provides such an approach. To sub- 
stantiate this claim, we will investigate two nontrivial 
problems: the second-order shifts in the critical temper- 
ature and chemical potential of the homogeneous Bose 
gas, and the condensation temperature of the Bose gas 
with delta-correlated disorder (22t l23| . 

The purpose of these studies is fourfold, (i) As already 
mentioned above, we would like to demonstrate that this 
RG method provides a transparent, elegant and efficient 
way to investigate critical thermodynamic properties of 
interacting Bose gases. First of all, quantitative results 
for, e.g., the critical temperature and chemical potential 
are readily obtained by numerically solving simple RG 
flow equations. The accuracy of the results is compara- 
ble to the other recent approaches mentioned above. Fur- 
thermore, the approach also allows one to obtain some 
qualitative insights into the essential behavior of criti- 
cal quantities. To illustrate this point, it will be shown 
how the nonanalytic dependence of the critical chemi- 
cal potential on the scattering length can be extracted 
analytically from the flow equations. 

(ii) We present the first second-order calculation of 
the critical temperature and chemical potential that 
does not rely on perturbative matching to the three- 
dimensional 0(2) scalar field theory. Most calcula- 
tions of c\ and, to the author's knowledge, all previ- 



ous evaluations of c^ 0, E3, IS 13 113 proceed by 
determining nonperturbative parameters of the classical 
three-dimensional (3D) theory and relating them to the 
quantum-mechanical problem through the matching pro- 
cedure of Refs. |j,|al24|- It would thus be desirable to 
evaluate the critical temperature shift in a way which 
is independent of this procedure. The present approach 
provides such an independent derivation since it works 
with the full 3 + 1-dimensional action and does not make 
use of dimensional reduction. 

(iii) By inv estig ating the Bose gas with delta-correlated 
disorder [22, |23, we show that the RG approach can 
easily be extended to the study of more complicated 
systems. This model problem has been of interest re- 
cently in order to understand qualitative features of cer- 
tain systems in random environments, such as He4 in 
porous media. The RG flow equations are derived with 
the help of the replica method. It is shown that in the 
high-temperature limit they reduce to the flow equations 
describingclassical spin systems with random exchange 
disorder Hi |H E3 Si El • 

(iv) Using the random fixed point of the flow equations, 
we calculate the condensation temperature in the pres- 
ence of particle interactions and disorder. The only deter- 
mination of T c so far [23J is of perturbative nature and 
neglects the atomic interactions. However, it is known 
that the interactions are an essential part of the model 
since they stabilize the disordered Bose gas against col- 
lapse 22]. They also play a crucial role in the critical 
flow obtained from the RG equations. It would thus be 
more satisfactory to calculate the condensation tempera- 
ture within an approach that directly includes the inter- 
actions. This is accomplished by our RG method. 

The approach presented in this paper is based on the 
application of the momentum-shell RG method [2l|] to 
a weakly interacting boson system. Since the 1970s, a 
number of studies have considered this question (see, 
e.g., Refs. [M EJ El El El El El) Whereas some 
of these works mainly focussed on the critical exponents 
|30l l3lL l35| . others used the RG equations to investi- 
gate critical thermodynamic quantities. In Ref. E3j the 
critical temperature of a two-dimensional (2D) Bose gas 
was examined, while Refs. J34| and J3|j studied critical 
properties of the 3D gas using the RG equations for the 
condensed phase. These equations, however, do not yield 
the correct behavior for, e.g., the critical temperature 
|37j . After clarifying the relation between the flow equa- 
tions for the condensed and the uncondensed phase |38| . 
the latter were successfully used to investigate the critical 
temperature of trapped Bose gases |39|, [4(J . 

The particular appeal of the RG approach lies in its 
conceptual simplicity combined with the fact that the 
thermodynamic properties of interest are easily extracted 
numerically from the flow equations. As discussed in de- 
tail in Refs. [3!j and |ij, the approach has some short- 
comings: it is only approximate, and some of the under- 
lying assumptions and approximations are not fully con- 
trollable and difficult to improve systematically. Never- 



theless, these problems are outweighed by the advantages 
of the method. 

We first consider the ordered homogeneous Bose gas. 
Its partition function can be written as a functional in- 
tegral 



Z= I V<t> I V<j>e 



is[<p,r] 



(2) 



with the Euclidean action 



S[<M* 



J dr /"d 3 xJ0*(r,x) 



or 2m 



-/' 



r,x) + ||0(r,x)| 4 } , 



(3) 



/i denoting the chemical potential, (5 the inverse tempera- 
ture, m the particle mass, and g the coupling constant for 
the s-wave contact interaction. The bosonic fields 0(x, r) 
are periodic in r with period H/3 and obey spatial periodic 
boundary conditions. We Fourier-decompose the bosonic 
fields according to 



(x ' r) = (hpvy/ 2 £ ak -" exp[i(k ' x 



-)], (4) 



with u> n — 2irn/hf3 and the cutoff condition |k| < A, 
and split them into a slow and a fast field component 
</>(x,t) = < (x, t)+i/) > (x, r), where the fast field <f> > only 
contains momenta within a thin shell A — dA < |k| < A 
near the cutoff. The Euclidean action is then expanded 
up to second order in </> > . Carrying out the functional 
integration over </> > , expanding the result up to second 
order in the coupling strength, and using lowest-order 
derivative expansion, one can write the partition func- 
tion in the original form (J2J, 10, where the integration 
is now over the slow field only and the parameters // and 
g have received an (infinitesimal) renormalization. Note 
that this approach retains the contribution of all Mat- 
subara frequencies and thus the full quantum- mechanical 
character of the problem. After rescaling all quantities to 
restore the original cutoff, the procedure is repeated un- 
til the bosonic fields are integrated out completely. The 
concomitant change of [i and g is described by the flow 
equations 



dM 

~dT 

dG 

~dl 



= 2M 



-7T = G- 
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2tt 2 



AbN BE {N B E + l) + 



1 + 2N BE 



2(E> - M) 



(5) 



(0) 



for the scaled quantities M = /3aMi G = A 3 f3\g/(3, 
b = 0//3a with the inverse cutoff temperature (3\ = 
m/h 2 A 2 , the scaled Bose-Einstein distribution Nbe — 
{exp[6(-E> — M)] — l}^ 1 and the scaled cutoff energy 
E > = 1/2. The flow parameter I describes the reduction 
of the original cutoff according to Ae~ l . Details of the 
derivation of Eqs. and © and further discussions can 
be found in Refs. |34l I3^l39ll4l| . Similar equations were 
also used, e.g., in Refs. I30l l3ll l32l 133 



The flow equations JSJ and © possess a nontrivial 
critical point at (M*,G*) = (1/12, 5tt 2 /72). Flow tra- 
jectories asymptotically reaching the fixed point repre- 
sent a system at the phase transition, and their initial 
values [M(0),G(0),b(0)] correspond to critical thermo- 
dynamic parameters. The initial value G(0) is related 
to the s-wave scattering length via a A = G(0)b(0)/[47r + 

2G(0)6(0)/tt] 0111. 

Using this framework, we now investigate the criti- 
cal temperature in the limit of small a = a/Ay with 
At = (2nh 2 P/m) 1 / 2 the thermal wavelength. Instead 
of calculating T c directly, we initially work out the shift 
8 = (nAy) c — C(3/2) in the scaled critical density. This 
quantity is determined by first of all calculating the ther- 
modynamic potential 



Who 



2^ 2 6(0) y 



(11 



-31 



In 



_ e -Ul)[E > -M(l)] 



(7) 

for a critical trajectory, taking the derivative Shom = 
— dWhom/dM(0)\ b ,o) g(o) an( i finally evaluating 8 — 



[2tt6(0)] 



3/2 



•She 



-C(3/2) |22| . The behavior of the critical 



temperature is then obtained through an inversion pro- 
cedure U . In addition, we will examine another nontriv- 
ial critical quantity, i.e., the (scaled) chemical potential 
v = (J3u) c = M(0)6(0) which was previously studied in 
Ref. (H. 

It is interesting to note that in Ref. J231, the (super- 
fluid) transition temperature of a two-dimensional inter- 
acting Bose gas was calculated by a related RG method. 
This calculation was also based on the evaluation of the 
critical density as the derivative of the free energy. Some 
analytic arguments using flow equations similar to Eqs. 
(|3|) and 10 were applied to work out the derivative and 
in this way to finally arrive at a qualitative estimate for 
T c . In the present case, the flow equations JSJ) and © 
can also be used to obtain some qualitative insights into 
the essential behavior of critical quantities. As discussed 
in the Appendix, one can easily show that they imply a 
functional form 



v(a) 



b[ v) a + b ( 2 v) a 2 \na + bi v) a 2 



for the dependence of the critical chemical potential on 
the scattering length. This is in agreement with the result 
of Ref. 0. 

Our ultimate goal, however, is to perform a quanti- 
tative calculation of critical properties. Determining 5 
and v analytically through Eqs. ©~0 is rather diffi- 
cult since one would have to very accurately relate initial 
conditions to the asymptotic behavior across a highly 
nonlinear regime of propagation. We thus perform our 
studies via a numerical analysis. Careful numerical eval- 
uation allows us to determine 6(a) and v(a) for values 
of S less than 10 -6 . The data obtained in this way is 
then used to determine the functional behavior of 8 and 
v up to second order in S through a fitting procedure. 
The fitting is somewhat complicated by the fact that the 



individual numerical errors for 8 and v do not appear 
to be independently Gaussian distributed so that a stan- 
dard x 2 analysis cannot be performed. We thus adopt 
the procedure described in the following. 

We assume a certain functional form for v and <5, e.g., 



8(a) = 


= b[ 5) a + b { 2 S) a 2 lna + b { 3 S) a 2 , 


(8a) 


v(a) = 


-- b[ v) a + b 2 v) a 2 \na + b { 3 v) a 2 , 


(8b) 



and perform a least-square-fit using all available data 
points below an upper limit So. Carrying out the fit 
for different values of So, we obtain coefficient functions 
b\ (So). This approach allows us to deal with a number 
of issues. 

(i) The actual functional form of v(a) and 8(a) also 
contains terms of higher order in a that we are neglecting. 
For this reason, the functions bi(ao) become dependent 
on &Q. Taking the limit So — > 0, however, the higher- 
order terms get less important, and we can read off the 
final values for bi as lima ^o bi(do). Furthermore, we 
find that the inclusion of higher-order terms in the fitting 
scheme can significantly reduce the So dependence of the 
functions bi, whereas the limiting values hardly change. 
However, we did not attempt to systemically investigate 
such higher-order terms. 

(ii) The numerical data cannot be obtained with arbi- 
trary precision but contain errors. This puts limitations 
on the practical feasibility of the procedure described 
in (i). For v(a) this issue is not severe, since the ini- 
tial conditions M(0) and b(0) can easily be calculated 
to very high accuracy by means of a bisection routine 
|39j . For <5(a), however, one has to evaluate the deriva- 
tive of the thermodynamic potential J7J which adds a 
significant amount of numerical noise. The noise shows 
up in an erratic behavior of the functions b\ (So) for 
small arguments. Only for larger So, the functions be- 
come sufficiently smooth. We read off the limiting values 
linia ^o fri(So) horn the smooth part, but have to allow 
for significantly larger errors bars of the coefficients bi for 
8(a) as compared to v(a), especially for the second-order 
terms. To some extent, the problem can be alleviated by 
using very large data sets. Error estimates for the coef- 
ficients b\ are obtained from inspecting graphs of the 
functions &j(do). 

(iii) The functional form of v(a) and <5(S) cannot be 
assumed to be known a priori [although the results of 
the Appendix suggest relations like ©], and from the 
numerical analysis alone, one is not able to prove unam- 
biguously that a certain form is indeed the correct one. 
Nevertheless, a useful fitting procedure should help to 
distinguish between more and less plausible behaviors. 
To show that our method is indeed able to this, we note 
that, for example, a slight change in the fit function from 
S 2 lnS to a 1 ' 99 In a has a strongly adverse effect on the 
behavior of the coefficients of (So) for v, i.e., they no 
longer smoothly reach a limiting value for So — > 0. By 
performing similar tests with varying functional forms, 



we have verified that the behavior indicated above in- 
deed provides a very good representation of the data. 
Finally, we also mention that the fitting procedure was 
successfully tested with sets of simulated noisy data. 

The results of the numerical analysis can be summa- 
rized as follows. We find that the behavior of v(a) and 
5(a) for small a, as obtained from our RG approach, is 
well compatible with a functional form of the kind JSJ, 
whereas other suitable forms could not be identified. For 
the coefficients wc find the values 



with 



b\ 



(s) 



°2 

u(S) 

°3 

b (v) 
°2 
b (v) 
°3 



-5.303 (±0.005) [-7.1], 

-177 (±5) [-146.81], 

-720 (±20) [-587], 

10.4495 (±0.0003) [10.4495], 

122.2 (±0.2) [100.53], 

425.5 (±0.7) [386.71], 



with the numbers in square brackets giving the values of 
Refs. [i( and [24(. Using the inversion formulas of Ref. 
[H , the coefficients of Eq. JTJ) are obtained as 



d = 0.9826 (±0.001) [1.32], 



= 23.8 
= 91 



(±0.7) 
(±3) 



[19.75], 

[75.7]. 



We find that the deviations of RG values from the re- 
sults of Refs. and range from 10% to 26% with 

the largest differences occurring for the coefficients b x 
and c\ . The accuracy of these results is thus comparable 
to some of the other recent approaches mentioned above. 
Earlier approaches even had deviations of far more than 
100% from the current best values [lj. Furthermore, it 
should be emphasized again that our calculation of the 
second-order shifts is completely independent from pre- 
vious works such as Ref. [J and does not make use of 
dimensional reduction. It also appears to be consider- 
ably simpler than these approaches. 

We now turn to the calculation of the condensa- 
tion temperature in the presence of disorder and as- 
sume the Bose gas to be exposed to a random poten- 
tial V(x) with average (V(x)) — and delta correlations 
(V(xi)V(x2)) — Rq5(x 1 — x 2 ). This model was stud- 
ied, e.g., in Refs. [22| and [23j to describe He4 in porous 
media. The parameter Rq is a measure for the strength 
of the disorder. To evaluate the partition function of the 
Bose gas we apply the replica method and consider an en- 
semble of N copies of the system. The partition function 
for this ensemble, averaged over the random potential, is 
given by plli^ 



ZN ) = {IT / V ^ j Wc 



-S W [*',*]/6 



(9) 



S W [0*,0] 



m> 



(It 



N 

Q = l 



U X ! T ) 



H- A 

dr 2M 



M 



^(x,r) + |C 2 (x,r)^(x,r) 



dr / dr' 
o Jo 

N N 



Rq 

2h 



d D o 



Yl J2 C( X : r )^a( X > r )C'( X , T ')<M X ,' 



(10) 



a— 1 a' — l 



Again making use of the first-order derivative expansion, 
the RG approach outlined above can be extended to deal 
with the additional terms in Eq. (|1(J|) introduced by the 
random potential. The resulting flow equations now also 
describe the renormalization of Rq. They read 



dM 

~~dT 



dG 

~dT 



dR 

~df 



2M 



1 
~ 2^2 
R 



-2GbN BE + NbRN BE 



1/2- M 
G + -^ [~G 2 b [AbN BE {N BE 
GR 



1) 



1 + 2N BE 



R 



2(1/2- M) 
1 

+ 2^ 

2 



6 



-R 



(1/2 -M) 2 
(NbR - 4bG)RbN BE {N BE 
4 



(11) 



(12) 



1) 



(1/2 -Mf 



(13) 



with the scaled disorder parameter R = A 3 (3\R n . Note 
that the number N of copies of the Bose gas now enters 
as a parameter into the flow equations. Within our one- 
loop approximation, the averaged partition function for 
the ensemble is given by 



(Z N ) = J] (1 - exp{-/?[ £k - /i(|k|; N)]}) 



-N 



(14) 



with e k = h 2 k 2 /2m and /x(|k|; N) = e~ 2l M(l; N)/f3 A the 
renormalized chemical potential with the trivial scaling 
removed. The parameters k — |k| and I are related by 
k = Ae~ l . The scaled free energy W = (/3a/ (3) (In Z) of a 
single copy of the system is obtained through the replica 
limit 



W = & hm ^ 



1 



[3 N^O 

1 



2tt 2 6(0) 



N 



die' 61 In [1 - e 



(15) 

-b(l)[l/2-M(l;N=0)h 



In other words, the free energy is calculated by integrat- 
ing the flow equations at N = [compare with Eq. (JJ]. 
Analogous types of disorder, i.e., the so-called random 
exchange or random temperature defects, have been thor- 
oughly examined in the context of classical spin systems 
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FIG. 1: Critical trajectories of flow equations l|110 - l|13fl pro- 
jected onto M-G plane. The trajectories have a — 4.9 x 1CP 4 
and g = 2.2x 10~ 4 (dashed-dotted), 4.0 x 10" 2 (dashed), 0.44 
(full), 1.78 (dotted), 20.0 (inset). 
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FIG. 2: Shift of critical temperature as function of g for fixed 
a = 4.9 x 10~ 4 . Full curve: RG result; dashed: Eq. ITrJl . 
Insets: (a) RG result after removing offset at g — 0; (b) shift 
of T c for o = l.l X 10" 2 . 



[H |2l HI HI ig. In the Ginzburg-Landau model, 
these defects are described by adding a random term 
<5r(x)S 2 (x) to the Hamiltonian, where S'(x) is the spin 
distribution and <5r(x) is delta correlated just like the 
potentials V(x) above 29]. It is therefore not surpris- 
ing that in the high-temperature limit b(l) — » the RG 
equations l|ll|l -l|13 (l reduce to the flow equations for the 
classical system (compare with Eqs. (2.9) of Ref. J2{|). In 
particular, the fixed-point structure of Eqs. (|1 1|) (|13|) is 
the same as for the spin system. For N = 0, these equa- 
tions possess three nontrivial asymptotic fixed points: 

1. (M?,Gl,Rl) = (1/12,5tt 2 /72,0) with eigenvalues 
((3 + \/249)/10, (3 - \/249)/10, §). This is the Wilson- 
Fisher fixed point well known from the ordered system. 
However, due to the presence of disorder it becomes 
destabilized and no longer describes a true phase transi- 
tion. This can be seen from the fact that two eigenvalues 
are positive. 

2. (M 2 *, G|, #|) = (3/22, 16tt 2 /121, 8tt 2 /121) with eigen- 
values ((1 + y/V77)/8M. - v / 177)/8, -i). This so-called 
random fixed point 25, (2y| is associated with condensa- 
tion in the presence of disorder. 

3. (M£,Gl,R%) = (1/18,0, -2vr 2 /9). This fixed point 
does not have a physical meaning: it has G$ — 0, but 
there is no phase transition in the disordered interaction- 
free Bose gas. 

FigureniUustrates the influence of disorder on the flow 
described by Eqs. (|llf) - H13[) . and shows the crossover be- 
tween the fixed points (FPs) 1 and 2. The figure dis- 
plays critical trajectories at fixed a — a/ At — 4.9 x 10~ 4 
and increasing g = (m/3/iiTh 2 )(Ro/a) = 6(0)i?(0)/47raA. 
The parameter g is a convenient dimensionlcss measure 
of the strength of disorder relative to the particle inter- 
actions. For each value of g, two very close numerical 
approximations i and i', i = 1, . ..,4, to the actual crit- 
ical trajectory are shown that eventually depart to the 



left and the right, respectively. 

For g -c 1 (trajectory set 1 in Fig. QJ, the critical 
flow is still essentially dominated by FP 1. Although 
this fixed point is no longer asymptotically attractive, 
the influence of disorder is so weak that we were not 
able to numerically find trajectories that go beyond FP 1 
and approach FP 2 (although we expect such trajectories 
to exist: if a critical trajectory remains in the vicinity 
of FP 1 long enough, the influence of disorder should 
eventually drive it away from FP1 and towards FP2). 
For increasing disorder (trajectory set 2 in Fig. 0, the 
critical trajectories go beyond FP 1, but we still cannot 
find trajectories reaching FP 2. This is only achieved 
for g m 1 (see set 3 and 4). In the regime of g ^> 1, 
the flow trajectories display "interaction-free" behavior 
for a long integration time (see inset). Only when M 
has almost reached 0.5, the influence of the interactions 
becomes noticeable and bends the trajectories towards 
FP 2. In this case M(l), G(l) and R(l) become very 
large. Since our RG approach is based on expanding 
the Euclidean action in G(l) and R(l), it is probably no 
longer valid in this regime. We note that a very similar 
overall behavior of the flow is also obtained for larger 
values of a. From these results, we conclude that our RG 
calculation is applicable in the regime of £><1. 

The critical flow is characterized by the vanishing 
of the physical renormalized chemical potential [i.e., 
M(l)e — ► for I — > oo]. Fixed point 2 can thus be as- 
sociated with a crossover transition between the normal 
and a locally condensed phase [23. Using perturbation 
theory, in Ref. |23j the corresponding shift in the con- 
densation temperature due to disorder was estimated as 



T r -T, 



(0) 



T, 



(0) 



Rr, 



T c (0) m 3 
67r 2 n 



-ga 



8tt 



Tr. 



3C(3/2) V Tr (0) 



(16) 



with n the Bose gas density and Tc the condensation 
temperature of the ideal gas. In the perturbative calcu- 
lation, effects of particle interactions arc neglected [note 
that the product gd on the right-hand side of Eq. l|ltj|) 
is independent of a]. However, the study of the fixed 
points and the critical flow in the RG approach suggest 
an essential role of the particle interactions even in the 
presence of disorder. It would thus be more satisfactory 
to calculate the shift within a nonperturbative approach 
that takes the interactions fully into account. With our 
RG method, such a calculation can be performed in an 
efficient manner. 

In Fig. |21 we compare the RG calculation of T c , which 
proceeds analogously to the case without disorder, to the 
perturbative result ^^. The relative shift AT C /T° is 
shown at fixed a as a function of the strength of disorder 
g. We focus on the regime g<l where the RG approach 
is applicable. For g — > 0, the RG shift smoothly tends 
to the (finite) value obtained in the calculation for the 
ordered interacting gas. This indicates that, in contrast 
to the perturbative calculation, our approach is able to 
simultaneously deal with the concurrent effects of inter- 
actions and disorder. 

As we show in inset (a) of Fig. [3 if the T c shift at 
g = is subtracted from the RG result, both curves are in 
relatively good agreement with each other. In particular, 
we find that for small, fixed a, the shift of T c in the RG 
behaves as d\ g+0(g 2 ). The numerical results suggest 
that the slope d\ ' depends on the scattering length a, 

and for small a the ratio to the corresponding value af @T 
of the perturbative result l|16() can be approximated as 



APPENDIX 

In the following, the qualitative behavior of the critical 
chemical potential v = (f3fJ,) c in the RG description of the 
homogeneous Bose gas is worked out. It appears quite 
difficult to do this directly from Eqs. J5J and ©. How- 
ever, these equations can be simplified in a way which, 
on the one hand, does not modify the essential behav- 
ior of v but, on the other hand, allows to investigate the 
dependence on a = a/Xx analytically. In a first step, 
we use that typically M <C E > with E> — 1/2 the 
scaled cutoff energy, and neglect the M dependence of 
the Bose-Einstein distribution N BB and of the denomi- 
nator 2(E > —M) which appears in Eq. JJjJ. Furthermore, 
we disregard the nontrivial renormalization of G at high 
energies due to two-body scattering and thus drop the 
term —G 2 b/2E > from Eq. ©. The relation between the 
initial condition for G(l) and the scattering length a is 
then given by G(0) = 8n 2 a/ -\/27r&(0). The flow equations 
now read 



dM 

~dT 

dG 

dl 



= 2M 



= G- 



2^bN BE , 



G 2 b 



2tt 2 



AbN BE (N BE + 1) + 



Nbl 



(17) 
(18) 



with Nbe = [exp(6(Z)£>) - l]" 1 
(|18fl is simplified further by setting 



b(l)N BE = 




1/E > 



E> 
In a second step, Eq 

b(l) > 1, 
6(0 < 1. 



This modification, which has only a small effect on the 
flow of G, turns Eq. (jTH|l into 



4 RG > 



4' 



pert) 



0.7 -90a. 



We thus find that in the RG calculation, T c decreases 
slower with g than predicted by perturbation theory. For 
increasing a, the discrepancy becomes larger, as can be 
seen from inset (b) in Fig. 

In summary, we have applied the momentum-shell RG 
method to examine two nontrivial questions regarding 
the condensation temperature of interacting Bose gases. 
First of all, we have calculated the shift in T c (and the 
critical chemical potential) for the homogeneous Bose 
gas without disorder up to second order in the scatter- 
ing length. Secondly, we have presented a nonperturba- 
tive calculation of the condensation temperature of dis- 
ordered Bose gases thereby taking the influence of inter- 
actions properly into account. Our results show that the 
RG method is an efficient and versatile tool to study crit- 
ical properties of interacting Bose gases and presents an 
alternative to the approaches used so far. 

Stimulating discussions with B. Kastening, A. Pelster, 
G. Metikas, and G. Alber are gratefully acknowledged. 




(19) 



where q = 20/2-7T 2 . The solution to this equation is given 
by 



G(0 = 



G(0)e l 
G(0)e l 



l + oG(0)[e' 




(20) 



The critical initial condition M c (0) for the chemical po- 
tential is determined by the condition e~ 2l M(l — > oo) = 
0. From Eq. (U7J thus follows that 



MM = —, / 2G(l)b(l)N BE e- 2L dl 

i r°° ~ 
= 7TT / 2G(0)e l b{l)N BE e- 2l dl 
^ Jo 

( G(0)e l 



+ - 



2tt 2 Jiny/m V 1 + 9 G (°)( e/ 



G(0)e' )b(l)N BE e- 2l dl 



h + h- 



(21) 



The first integral I\ equals 6(0) 1 4('(3/2)a, whereas I2 
can be rewritten as 



gG(0) s 

2tt 2 



dx — 

1 



y/x-\ 



1 



g G(o)v^(oy(i-i/Vi)e^>-r 

(22) 
Expanding the second factor in the integrand as 
l/(e xE > - 1) = 1/xEy - 1/2 + • • • , one sees that ex- 
pression 122fl can be written as 



h = 



1 

W) 



—a 2 \na + b'^a 2 + 0(a 3 ) 



(23) 



It is interesting to note that, similarly to the effective- 
field theory calculation of Ref. [24J, only the coefficient 
of a 2 In a can be determined by the perturbation expan- 
sion. All other terms apparently receive contributions 
from all expansion orders and thus cannot be calculated 



in this way. Putting together the above results and using 
numerical data to obtain the coefficient of the a 2 term, 
we finally find 



640 

v = (/?/i) c = 4C(3/2)fi+ a 2 lna + 908a 2 + O(a 3 ) (24) 

■K 



in the simplified approach studied here. It should be re- 
marked that even in this simple model, the term linear 
in a has the correct prefactor, whereas the coefficients 
of the two second-order terms differ by only a factor of 
about 2 from the values given in ref. |24| and the results 
of our original RG equations. We also note that from 
the numerical simulation of Eqs. I|17l) and H19|) we can 
accurately reproduce the coefficient 640/7T of the a 2 In a 
term, thus confirming the accuracy of the fitting proce- 
dure described in the main text. 
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